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20.  ABSTRACT  (Continued): 

'It  was  found  that  modifications  to  the  usual  Morison  approach  are  required  to  adequately  ac¬ 
count  for  the  orbital  motions  of  the  fluid  and  to  account  for  the  orientation  of  the  orbits  with  re¬ 
spect  to  the  cylinder  axis.  The  consequences  of  these  findings  are  discussed  for  inclined  cylinders  in 
waves  and  for  cylinders  in  short-crested  seas.  The  axial  variations  of  the  wave  force  on  vertical  cyl¬ 
inders  are  examined  in  order  to  establish  error  bounds  for  the  common  practice  of  assuming  constant 
values  of  C^  and  CD  over  the  entire  span.  Also,  the  methods  of  computing  force  transfer  coeffi¬ 
cients  from  a  force  record  are  examined  and  errors  of  method  are  identified  and  briefly  discussed.  . 
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THE  SENSITIVITY  OF  WAVE  FORCE  COMPUTATIONS 
TO  COMMON  ERRORS,  UNCERTAINTIES, 

AND  HYDRODYNAMIC  APPROXIMATIONS 


1.0  INTRODUCTION 

In  the  face  of  rising  construction  costs,  the  increased  importance  of  dynamic  loadings,  and 
in  view  of  the  more  hostile  sites  under  consideration,  the  accuracy  of  wave  force  computations 
becomes  a  critical  question.  The  least  well  understood  wave  loading  regime,  hence  the  one 
with  the  least  accurate  descriptions,  is  the  regime  wherein  both  drag  and  inertia  forces  are 
important.  The  basis  for  most  wave  force  computations  in  this  regime  is  the  so-called  Morison 
equation.  In  the  past  the  Morison  approach  has  been  tailored  with  some  success  to  particular 
applications.  However  attempts  to  generalize  the  approach  have  not  been  successful  and  as  a 
consequence,  large  uncertainties  can  accompany  a  new  application.  The  discrepancies  between 
prediction  and  observation  are  often  as  large  as  50  to  100  percent  [1). 

For  the  most  part  these  inaccuracies  stem  from  a  poor  understanding  of  the  unsteady  vor¬ 
tex  flows  which  occur  in  this  regime  and  from  several  simplifications  inherent  in  the  usual  Mor¬ 
ison  approach.  In  lieu  of  a  complete  hydrodynamic  description  of  the  complex  vortex  wake 
flows,  an  unlikely  accomplishment  in  the  near  future,  an  improved  approach  to  wave  force 
prediction  appears  to  be  one  that  can  account  for  the  major  differences  between  wake  flows 
which  may  occur  for  various  wave-cylinder  combinations.  This  report  outlines  such  an 
approach  and  at  the  same  time  identifies  uncertainties  that  can  arise  from  differences  in  the 
methods  of  force  or  force  coefficient  calculation. 
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Most  investigations  of  the  fluid  forces  on  an  object  immersed  in  an  unsteady  flow  begin 
by  resolving  the  force  into  drag  and  inertia  components  [2].  It  is  less  widely  known  that  the 
general  unsteady  force  also  contains  a  history  term  and  should  be  written  as  the  sum  of  three 
contributions  in  the  form, 

frow/O)  -  F^u)  +  FlMftie(u)  +  FHlstwy(u.  lif.  r)  (I  D 

where  u  and  it  are  the  velocity  and  acceleration,  respectively,  which  depend  only  on  the  time  t. 

Some  years  ago,  Morison  et  al  [3]  made  use  of  this  approach  in  attempting  to  describe  the  wave 

forces  on  a  cylinder.  In  essence,  they  took  an  elemental  force  on  the  cylinder  to  be  of  the  form 

dF{t,s)  «  dFom  +  dFlnertla  (1.2) 

where 

dFo^-Mlp  D  CDu\u\ds  (1.3) 

and 

dFlntnla  "  »/4  p  D2  C„  it  ds  (1.4) 

and  with  s  measured  along  the  cylinder  axis.  Since  this  equation  first  appeared  it  has  become 

the  basis  for  the  majority  of  practical  wave  force  computations  and  is  widely  known  as  the  Mor¬ 
ison  equation.  The  unknown  force  transfer  coefficients  CM  and  C0  are  usually  taken  as  time 
invariants  so  that  the  unsteadiness  of  the  force  resides  entirely  in  the  variations  of  u  and  u. 

Along  with  the  recent  growth  of  offshore  construction,  considerable  effort  has  been 
expended  to  collect  values  of  Cu  and  CD  and  to  correlate  the  changes  in  the  coefficients  with 
the  important  parameters  of  the  problem.  Except  for  a  few  simple  cases,  notably  not  including 
progressive  wave  flows,  the  variations  in  the  force  transfer  coefficients  have  not  been  com¬ 
pletely  accounted  for.  To  be  sure,  one  can  obtain  reasonably  good  force  estimates  but  even 
with  great  care  in  the  analysis  errors  of  20  to  30  percent  in  the  total  force  are  common  and 
errors  in  the  local  force  are  often,  as  great  as  100  percent  or  more  ll). 
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The  objective  of  the  present  investigation  is  to  improve  the  resolution  of  the  Morison 
approach  by  first  examining  the  sensitivity  of  the  computed  wave  forces  (or  computed 
coefficients)  to  the  common  simplifications  of  method  and  to  the  common  hydrodynamic 
approximations.  Then,  alternatives  or  modifications  that  can  improve  the  resolution  are  con¬ 
sidered.  These  are  compared  to  wave  force  data  obtained  in  a  laboratory  channel  and  are  com¬ 
pared  to  previous  results  for  simple  oscillatory  flows  about  cylinders  [4,  5,  6,  7,  8],  By  way  of 
example,  the  present  approach  can  be  introduced  by  considering  a  simplification  of  method  and 
a  hydrodynamic  approximation  that  already  have  occurred  in  passing  from  Eq.  1.1  to  Eq.  1.2. 

The  simplification  is  due  to  the  apparent  lack  of  a  history  term  in  Morison’s  equation  as 
compared  to  the  general  Eq.  1.1.  However,  if  the  flow  is  purely  periodic  (0  -  «/)  the  omission 
does  not  actually  occur  because  the  now  periodic  history  term  is  included  in  the  periodic  drag 
and  inertia  terms.  This  can  be  demonstrated  as  follows.  Let  the  drag  and  inertia  forces  be  writ¬ 
ten  in  the  usual  way 

~  CD*  u\u\ds  (1.5) 

and 

^  -C^iids.  (1.6) 

With  little  or  no  loss  in  generality  the  periodic  history  term  can  be  expanded  in  the  form 

dFHamy  -  (1/2  p  D  CJ>H  it  |v |  +  ir/4  p  D2C\,h  u)  ds  (1.7) 

becaure  both  u  and  u  will  be  composed  of  orthogonal  periodic  functions.  One  can  then  redefine 

coefficients  to  obtain 

dF(9)  -  (1/2  p  D  CDu\u\+n/4p  D1  CM  u)  ds  (1.8) 

where 

CD-CI>  +  Ch,  (1.9a) 

and 
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Cv-  Cfc+Cfc*,  (1.9b) 

represent  composite  force  transfer  coefficients.  The  significance  of  this  simplification  is  that  the 

composite  coefficients  obtained  from  either  of  two  separate  experiments  may  not  be  appropriate 

when  the  two  flows  are  superimposed.  Important  problems  in  this  category  are  the  forces  on  a 

structure  exposed  to  a  wave/current  system  or  one  experiencing  oscillatory  motion  as  well.  The 

history  term  will  not  be  considered  further  in  this  investiption;  here  the  coefficients  CM  and 

Co  represent  composite  values  for  nominally  purely  periodic  flows. 

The  hydrodynamic  approximation  in  Eq.  1.2  is  the  assumption  that  Eq.  1.1  can  be 
extended  to  two-dimensional  flows  by  taking  a  differential  approach  and  then  assuming  that  the 
distributed  force  depends  only  on  the  local  kinematics  of  the  flow.  In  fact,  the  variation  of  flow 
along  the  span  will  introduce  three-dimensional  effects  in  a  number  of  ways.  The  axial  pressure 
gradient  will  produce  an  axial  component  of  velocity  which  can  greatly  influence  wake  forma¬ 
tion.  Less  obvious  is  the  effect  arising  from  the  wake  which  may  be  swept  back  over  one 
segment  of  the  structure  after  being  generated  at  another  segment  under  different  flow  condi¬ 
tions.  Three-dimensional  effects  have  been  observed  to  have  a  pronounced  influence  on  the 
drag  in  steady  bluff  body  flows  [IS, 16].  This  approximation  is  examined  further  in  Section  3.0 
of  this  report  along  with  other  sources  of  errors  and  uncertainties.  Before  these  are  considered 
a  brief  review  of  previous  results  from  planar  or  one-dimensional  harmonic  flows  is  presented. 

2.0  SIMPLE  OSCILLATORY  FLOWS  -  A  REVIEW  OF  RESULTS  AND  METHODS 

The  experiments  of  Sarpkaya  [4,  S],  Bearman  et  al.  [6,  7],  Maull  and  Milliner  [8],  and  to 
a  slightly  lesser  extent  the  experiments  of  Keulegan  &  Carpenter  [9]  belong  in  the  category  of 
simple  oscillating  flows  over  a  cylinder.  Each  investigation  approaches  the  conditions  for  which 
Eq.  1.1  is  written  and  therefore  does  not  require  several  of  the  usual  assumptions  and 
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simplifications  of  Mori  son’s  equation  necessary  for  wave  flows.  Most  importantly,  each  flow  is 
one-dimensional  or  nearly  so. 

Sarpkaya’s  cylinders  were  mounted  in  the  horizontal  section  of  a  large,  U-tube  apparatus. 
A  similar  arrangement  was  used  by  Maull  and  Milliner  [81  and  Bearman  et  al.  (6,7],  Once  set 
in  motion,  the  fluid  oscillated  back  and  forth  over  the  cylinder  at  the  natural  period  of  the  water 
column.  Keulegan  &  Carpenter  placed  their  cylinders  horizontally  under  the  node  of  a  standing 
wave.  In  ail  of  these  studies  the  fluid  velocity  and  acceleration  fields  could  be  adequately 
described  by  the  single  components 

u  —  Um  cos®  (2.1a) 

and 

w  -  -  -y  Un  sin®.  (2.1b) 

In  each  case  the  flow  did  not  vary  along  the  cylinder  axis  and  Eq.  1.8  could  be  directly 
integrated  over  the  cylinder  length  L  to  obtain 

F(9)  -  1/2  p  DCdL  Vjt,  cos®  | cos 6 1  -  irJ/2 Tp  D 1  L  CM  UM  sin®  (2.2) 

for  the  total  unsteady  force.  By  matching  this  expression  with  the  experimental  results  it  has 

been  found  that  the  variations  in  Q>  and  CM  could  be  accounted  for  through  the  use  of  any 
two  of  the  following  three  parameters. 


UmT 

-  — . s  period  parameter  =  K  (Keulegan-Carpenter  number) 

(2.3a) 

UmD 

- =  Reynolds  numbers  =  R 

V 

(2.3b) 

D 2 

—  =  frequency  parameter  = 
vT 

(2.3c) 

Only  two  of  the  three  parameters  are  independent  because  R  -  p  ■  K.  The  period  parameter 
and  Reynolds  number  are  often  emphasized  because  of  the  analogy  with  inverse  Strouhal 
number  and  Reynolds  number  in  steady  flows.  Unfortunately  both  K  and  R  contain  Vm  and 
this  can  mask  important  frequency  effects.  The  use  of  the  frequency  parameter  p  has  a  certain 
practical  advantage  since  it  is  a  constant  for  a  given  wave/cylinder  combination. 
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Using  the  Keulegan-Carpenter  number,  the  above  equation  for  the  fluid  force  may  be 
written  in  the  non-dimensional  form 

/(#)  -  ~  Cd  cos«|cosfl  I  -  ~  CM  sin«,  (2.4) 

which  will  be  employed  in  the  balance  of  this  investigation.  In  this  form  it  is  also  readily 
apparent  that  the  ratio  of  the  peak  drag  force  to  the  peak  inertia  force  is  given  by  CDK/ir 2  CM 
and  this  should  be  kept  in  mind  when  assessing  the  relative  importance  of  errors  in  CM  and  CD 
for  various  ranges  of  K.  The  conclusions  that  have  been  drawn  from  the  one-dimensional  or 
planar-flow  experiments  are: 

i)  Morison’s  equation  appears  to  adequately  describe  the  in-line  force  for  simple  one¬ 
dimensional  oscillating  flows  about  circular  cylinders. 

ii)  Tabulations  of  CM  and  Q>  are  available  for  predicting  in-line  fluid  forces  on  circular 
cylinders  in  one-dimensional  oscillating  flows;  see  Figures  1  and  2. 

iii)  Several  shortcomings  in  the  Morison  approach  have  been  identified  for  non-circular, 
sharp-edged  cylinders  [6,  7], 

iv)  Large  cycle-to-cycle  variations  in  Cm  and  CD  can  occur  even  for  nearly  uniform  total 
forces  as  measured  by  the  root-mean-square  (rms)  value  [6,  7,  81. 

v)  A  side  or  lift  force  is  generated  which  may  be  nearly  as  large  as  the  in-line  force 
(9,101. 

As  stated  in  the  introduction,  the  direct  application  of  these  results  to  wave  force  predic¬ 
tion  has  resulted  in  large  errors.  This  is  not  surprising  in  view  of  the  added  complexities  of 
wave-induced  flows  but  it  suggests  certain  objectives.  For  example,  one  could  decide  to  retain 
the  usual  form  of  Morison’s  equation  and  then  catalog  values  of  CM  and  CD  for  each  different 
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situation.  This  is  the  pragmatic  approach  taken  by  much  of  the  offshore  industry  and  will,  or 
course,  pose  some  difficulty  for  every  new  situation.  Alternatively  one  could  modify  the  equa¬ 
tion  in  search  of  forms  that  permit  a  single  tabulation  or  several  tabulations  of  CM  and  CD  to 
certain  levels  of  accuracy  and  equation  complexity.  Both  approaches  are  examined  in  the  fol¬ 
lowing  sections,  but  first  it  is  instructive  to  present  the  three  common  methods  for  determining 
Cm  and  CD  from  a  force  record. 


The  first  method  makes  use  of  the  fact  that  the  drag  force  is  a  maximum  when  the  inertia 
force  is  zero  (®  ->  0,  n)  and  vice  versa  ( ®  —  ir/2,  in/2).  The  measured  force  is  then  examined 
at  these  special  times  in  the  cycle  and  set  equal  to  either  the  drag  or  inertia  force  as  appropriate. 
This  method  is  very  simple  and  has  the  added  advantage  of  relying  primarily  on  peak  forces. 
However,  it  lacks  averaging  to  cancel  spurious  errors,  and  in  more  complex  flows  the  presence 
of  harmonics  can  distort  the  results  (see  section  3.1).  Morison  et  al  13)  used  this  approach  but 
with  the  advent  of  digital  computers  to  handle  the  data  reduction  single-point  methods  have 
largely  beep  abandoned. 


The  second  method  employs  a  Fourier  analysis  of  the  force  coefficients.  Let  the  non- 
dimensional,  measured  force  be  given  by  Eq.  2.4,  and  then  multiply  it  by  sin®  and  integrate 
over  a  cycle.  This  yields  the  inertia  coefficient 


_  J„  /mwjW)  sin®  dO  K  r2. 

Cm  ~2  2tr  ™  T  J0  fmtm  sin®  d9. 

tX  ,m  * 


(2.5) 


When  the  measured  force  is  multiplied  by  cos®  and  integrated  over  a  cycle,  the  drag  coefficient 


is  given  by 


CD~ 


X  /mMJ(®)  cos 0  dO 

/► 

cos2®|cos®|d® 

0 


-i  rv 

1  Jo  J* 


cos®  dO 


(2.6) 
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The  integrals  can  be  evaluated  numerically  using  the  measured  force.  The  Fourier  approach 
emphasizes  the  periodic  nature  of  the  problem  and  can  be  extended  to  find  any  desired  har¬ 
monic  in  the  measured  force  record  as  was  done  originally  by  Keulegan  &  Carpenter. 


The  third  method  is  a  least-squares  curvefit1  given  by  the  weighted  sum  of  the  squares  of 
the  force  residuals, 


5- I  W—-/W 

i-i 


(2.7) 


subject  to  the  constraints 


bS  bS 


-  0. 


(2.8) 


BCU  bCo 

This  will  minimize  the  differences  between  the  measured  and  calculated  forces  in  the  least- 


square  sense.  Performing  the  differentiation  for  the  general  coefficient  C  gives 

H  -  Jc  I  WX fL»  -  2  /**,  m  +  -  0 


or 


I W, 


9/L 


bC 


-2  /, 


izm  +  2  m  iZM. 

meas  be  J  be 


-  0 


and 


(2.9) 


(2.10) 


,  W  ,  ftt  I  >/<>'> 

1  “ ac  1  ’•  f'>l>  SC 

where  the  summation  is  understood  to  be  over  the  N  data  points. 


(2.11) 


Introducing  Eq.  2.4,  to  get 


and 


d  f 

T7T-  “  COStfJcOStfjl 
vC  d 


a/, 

a  Or 


'See  Appendix  A  for  a  description  of  the  least-squares  method. 
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yields  two  equations 


I  W,  cos«*|cos0,| 
Cd“  I  W,  cos40,  + 


nJ  Z  Wj  cos9*|cosfl*|sinfl, 
~KC"  Z  W,  cos40, 


(2.14) 


and 


I  W,  fmtas ,  sin9*  l  wt  cos» ,|cos0,|sin0* 

Cu - - - L - +  CD - - - ! - ! - L  (2.15) 

—  £  W,  sin20,  ~  I  W,  sin20, 

K  A 

which  can  be  solved  for  CM  and  CD.  Clearly,  if  the  quantity  Z  W,  cose*  |cose*|sine*  is  zero 
then  the  solutions  for  Csi  and  CD  are  given  by  the  first  terms  of  Eqs.  2.14  and  2.15.  For  uni¬ 
form  statistical  weighting  ( W,  -  constant)  and  for  equal  increments  in  e*  from  0  to  2n 
inclusive,  this  quantity  is  identically  zero  and  the  solutions  for  CD  and  Cu  reduce  to 


and 


CD 


£  fmeas,  COS0,|cOS0,| 

i  cos4e* 


£  /m«»,  sine, 
n7  Z  sin20* 


(2.16) 


(2.17) 


Analogous  equations  can  be  obtained  using  an  integral  error  criterion2  with  the  results 


and 


(2.18) 


CM  -  — T  S  '  fmm  sine  <J0  (2.19) 

which  were  utilized  by  Sarpkaya  [5],  The  advantage  of  the  integral  formulation  is  that  the 
terms  containing  cos0)cos0lsin0  drop  out  during  the  integration,  but  it  must  be  employed  with 
some  caution  because  is  usually  a  discrete  function  and  the  integral  criterion  becomes 
strictly  inappropriate.  Under  the  conditions  set  forth  in  Appendix  A  the  integral  and  discrete 

least-squares  method  can  be  used  interchangeably.  Those  conditions  are  tantamount  to  the 
2See  Appendix  A 
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conditions  imposed  above  to  make  the  discrete  summations  of  cos0|cos0|sin0  vanish  and 
require  the  integrals  to  be  evaluated  in  a  certain  way  as  well. 

Upon  comparing  the  integral  least-squares  results  and  the  Fourier  results  the  expressions 
for  Cw  are  found  to  be  identical  but  the  two  methods  produce  different  equations  for  CD.  This 
difference  can  be  examined  by  first  supposing  that  the  measured  fluid  drag  force  is  proportional 
to  cos0|cos0|.  For  this  condition  the  two  methods  produce  the  same  numerical  result  as  they 
should.  On  the  other  hand  a  drag  force  proportional  to  cos0  or  proportional  to  cos30  will  pro¬ 
duce  about  a  4  percent  difference  between  computed  coefficients.  Sarpkaya  applied  both 
methods  to  his  measured  forces  and  found  a  small  but  consistent  difference  between  results 
which,  according  to  the  above,  would  imply  the  presence  of  large  even  harmonics  of  the  force 
other  than  those  due  to  cos|cos0|.  In  other  words,  Morison's  equation  may  not  account  for  all 
of  the  significant  in-line  forces  even  for  a  simple  planar  flow. 

If,  for  some  reason,  the  statistical  weighting  is  not  uniform  or  a  few  irregular  intervals  in 


0,  are  chosen,  then  the  full  Eqs.  2.14  and  2.15  must  be  solved.  To  do  this  we 

first  make  the 

following  substitutions 

f M  “  fmt as 

(2.20a) 

ft- Y  sin0i 

(2.20b) 

fp  -  COS0,|COS0,| 

(2.20c) 

in  Eqs.  2.14  and  2.15  to  get 

„  „  I  W,fDf, 

D  ZW,fh  M  I  W,fl 

(2.21) 

and 

„  -I  WfM,  .  „  I  W,fDf, 

M  I  w,  f}  D  1  W,  f}  ' 

(2.22) 
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The  solutions  for  the  coefficients  are 


I  W  f}  -  I  W,  f,  fD  1  Wt  f,  fM 
°  ’  'LWlf}  'LWfl-(ZWt  f,  /0)2 

and 


r  ~r  *  fl**  Wififp  n 

These  equations  are  of  the  form  given  by  Chakrabarti  et  al  (10,  11]  in  their  studies  of  wave 
forces  on  inclined  cylinders.  The  sign  changes  in  Eq.  2.24  as  compared  to  Chakrabarti's  expres¬ 
sion  are  a  consequence  of  the  minus  sign  before  the  inertia  term  in  Eq.  2.2  which  is  in  turn  a 
result  of  the  definitions  in  Eq.  2.1. 


One  important  source  of  error  which  is  common  to  all  of  the  above  methods  is  that  due  to 
an  uncertainty  in  the  phase  between  the  fluid  motion  cycle  and  the  measured  force  cycle.  This 
uncertainty  can  arise  from  a  phase  shift  in  some  portion  of  the  sensing  or  signal  processing  sys¬ 
tems  or  may  occur  physically  as  a  result  of  cycle-to-cycle  variations  in  the  wake  processes. 
Whatever  the  cause,  the  sensitivities  of  the  usual  Morison  coefficients  to  an  arbitrary  phase 
shift  are  given  by 

8w  C 

Cm! Cm  —  .  J*  sin  £  +  cos  $  (2.25) 

SK  C m 

and 

Cp/Cp  “  /(4)  _  ~  sin  4  (2.26) 

when  the  measured  force  is  assumed  to  be  given  by  Morison’s  equation  but  phase-shifted  in 
the  integral  least-square  expressions  by  an  amount  The  primes  in  the  numerators  are  used 
to  denote  the  phase-shifted  coefficients  as  compared  to  the  ’correct”  values  in  the  denominators 
and  on  the  right-hand  sides.  The  term  denoted  by  /(*)  is  a  lengthy  trigonometric  expression 
which  is  closely  approximated  by  cos  4  for  the  present  purposes.  Two  features,  which  are 
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expected,  become  apparent  and  quantified  in  theae  expressions.  First  is  the  tendency  to  swap 
one  component  for  the  other  as  the  phase  shifts  [11  and  the  second  is  the  marked  influence  of 

CD 

the  quantity  —  —  which  is  the  ratio  of  peak  drat  force  to  peak  inertia  force  as  discussed  in 
*  Cv 

connection  with  Eq.  2.4.  These  equations  are  plotted  in  Figure  3  for  several  values  of  the 
phase  shift  <fi.  In  viewing  the  figure  it  must  be  kept  in  mind  that  large  errors  in  either  C*  or 
CD  are  relatively  unimportant  when  the  other  force  component  dominates  the  total  force  as 
w1  Cd 

determined  by  the  ratio  —  — r~.  The  most  striking  results  in  the  figures  are  the  significant 

A  C*/ 

errors  C5-10  percent)  which  can  occur  for  phase  shifts  as  small  as  1  or  2  percent  of  the  cycle. 
Clearly,  special  care  must  be  taken  to  remove  or  identify  phase  shifts  in  the  measurement  sys¬ 
tem. 


Before  considering  wave  flows  one  further  comment  is  in  order.  Experiments  with  simple 
oscillatory  flows  as  described  here  have  revealed  the  presence  of  large  transverse  or  side  forces 
that  are  not  accounted  for  in  the  Morison  equation.  The  side  or  lift  force  as  it  is  sometimes 
called  arises  from  the  asymmetric  vortex  shedding  that  occurs  for  moderate  and  large 
Keulegan-Carpenter  numbers.  If  K  is  large  enough,  the  transverse  force  can  oscillate  a  number 
of  times  during  a  wave  cycle  in  relation  to  the  number  of  vortices  shed.  This  not  only  raises 
the  possibility  of  flow-induced  vibrations  but  also  raises  the  possibility  that  the  resultant  peak 
force  on  the  cylinder  will  be  greater  than  the  peak  force  predicted  by  Morison's  equation.  Max¬ 
imum  side  forces  of  comparable  magnitude  to  the  maximum  in-line  force  have  been  recorded 
in  one-dimensional  experiments.  This  difficulty  may  not  be  as  severe  in  wave  flows  because  the 
three-dimensional  character  of  such  flows  tends  to  limit  the  span  wise  coherence  of  the  shedding 
processes.  It  should  be  noted,  however,  that  cylinder  vibrations  can  restore  the  coherence  and 
therefore  deserve  consideration. 
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3.*  WAVE  FLOW  APPROXIMATIONS  AND  UNCERTAINTIES 

The  velocity  and  acceleration  fields  induced  by  a  progressive  surface  wave  are  nominally 
two-dimensional  with  components  (u,  w)  and  (u,  w),  respectively.  The  w  components  coincide 
with  the  direction  of  wave  advance  while  w  and  w  occur  in  the  vertical  direction.  Short-crested 
or  multi-directional  seas  will  produce  additional  components  of  velocity  and  acceleration.  In 
order  to  apply  Morison’s  equation  to  such  a  complex  flow  field  some  adaptations  and 
simplifications  are  required.  This  section  considers  the  common  approximations  and 
simplifications  in  order  to  assess  their  influence  on  the  accuracy  of  predicted  wave  forces  and  to 
assess  their  effects  on  computed  force  coefficients. 

The  additional  hydrodynamic  complexities  that  are  encountered  in  wave  Bows  may  be 
summarized  as  follows. 

e  The  flow  is  not  always  perpendicular  to  the  cylinder  axis. 

e  The  wake  is  not  always  swept  back  and  forth  over  the  cylinder. 

•  The  flow  is  not  always  uniform  along  the  span. 

o  The  flow  velocities  and  accelerations  are  not  explicitly  known. 

The  first  two  effects  are  largely  determined  by  the  eccentricity  of  the  water  particle  orbits  and 
by  the  orientation  of  the  cylinder  with  respect  to  the  orbits  [12].  The  variation  of  the  flow 
along  the  span  will  introduce  three-dimensional  effects  in  a  number  of  ways.  Obviously,  the 
instantaneous  velocities  and  accelerations  will  have  an  axial  variation  which  can  alter  the  flow 
forces  away  from  the  distribution  predicted  using  one-dimensional  results,  loss  obvious  is  the 
three-  dimensional  effect  arising  from  the  wake  which  may  be  swept  back  over  or  near  one  seg¬ 
ment  of  the  cylinder  after  being  generated  at  another  segment  under  different  flow  conditions. 
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Three-dimensional  effects  have  been  observed  to  have  a  pronounced  influence  on  the  drag  in 
steady  bluff  body  flows  (13,  14]. 

The  flow  velocities  and  accelerations  are  usually  computed  from  a  particular  wave  theory 
which  has  been  matched  to  the  wave  shape.  Since  different  wave  theories  will,  in  general,  pro¬ 
duce  different  pairs  of  force  transfer  coefficients  for  the  same  data,  it  is  often  suggested  that  the 
user  be  consistent  in  the  selection  of  wave  theory  [1].  Many  investigators  have  employed 
higher-order  expenaions,  e.g.  —  Stokes  Sth  order,  in  wave  force  computations  and  report 
improved  matching  to  the  wave  force  [2].  In  all  likelihood  these  improvements  are  largely  a 
result  of  the  additional  terms  in  the  curve  fitting  function  and  bear  little  relation  to  any 
improvements  in  the  hydrodynamic  description  as  discussed  in  the  following.  In  many 
instances  the  usual  approach  reduces  to  a  selection  of  a  wave  theory  that  best  matches  the 
observed  wave  surface  profile  or,  in  other  words,  the  velocity  and  acceleration  fields  are 
obtained  by  matching  computed  and  observed  displacements  at  the  free  surface.  Since  the  velo¬ 
city  and  acceleration  are  essentially  derivatives  of  displacement  it  is  clear  that  a  good  match  to 
the  surface  displacements  does  not  insure  an  equally  good  match  to  the  interior  kinematics  of 
the  flow.  Stokes  expansions  are  routinely  invoked  for  describing  the  kinematics  of  finite-height 
waves,  so  that  it  is  constructive  to  examine  their  use. 

First,  the  dynamical  surface  boundary  condition  can  be  used  to  illustrate  the  problems  of 
obtaining  velocities  from  a  matching  of  surface  displacements.  Let  the  surface  elevation  v  be 
expanded  in  the  regular  perturbation  series 

V  “  +  m«2  +  nj*2  +  0(«4)  (3.0 

where  «  «  1  is  the  wave  slope.  The  velocities  and  accelerations  can  be  similarly  expanded  but 

not  to  the  same  order  of  accuracy  if  the  dynamical  boundary  condition  given  by 
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it  to  be  uniformly  satisfied.  In  particular,  the  quadratic  terms  in  the  velocity  require  that  a 
remainder  of  order «"  in  the  surface  displacement  will  correspond  approximately  to  a  remainder 
of  order  «"/}  in  the  velocity  expansion.  Hence,  the  use  of  such  expansions  will  require  a 
disproportionately  large  number  of  terms  to  significantly  improve  the  representations  of  velocity 
and  acceleration  beyond  their  linear  descriptions.  In  view  of  the  ease  of  superposition  of  linear 
waves  and  the  practical  importance  of  directionality  and  short-crested  seas,  higher  order  expan¬ 
sions  or  theories  hardly  seem  worthwhile  except  possibly  in  the  extreme,  single  design  wave 
approach.  The  use  of  Stokes-type  expansions  can  also  introduce  other  sources  of  error  and 
uncertainty  as  described  below. 


The  two  problems  of  interest  are  the  use  of  Morison’s  equation  to  predict  wave  forces  and 
the  use  of  N^rison's  equation  to  obtain  the  force  coefficients  Cv  and  CD  from  a  wave  force 
record.  To  begin,  let  the  velocity  expansion  be  written  symbolically  as 

u  “  £  Un  cos  n».  (3.3) 

H-l 

so  that  the  acceleration  is  given  by 

“  *  Hr  I  U.  sin  lid.  (3.4) 

1  »-i 

where  r  represent  the  order  of  the  expansion  and  •  «  tor  -  wt.  This  notation  is  adopted  for 
convenience  and  it  should  be  noted  that  the  forms  of  the  coefficients  UK  depend  on  both  n  and 
r.  To  simplify  the  discussion,  the  usual  Morison  approach  is  taken  for  a  horizontal  cylinder  in 
waves  so  that  so  the  nondimensional  force  becomes 


<■  '  U.  U  i 

/,<•) -c#n  —in 

«-l  J-\  t M 


cos  n  9\co»J  •!  +  n|— -|  **n  "• 


A— I 


(3.5) 


when  the  characteristic  velocity  Um  is  equal  to  U i,  a  constant.  As  expected,  a  large  number  of 
harmonics  are  rapidly  introduced.  For  fixed  values  of  the  coefficients  C*  and  CD  the  presence 
of  these  harmonics  can  be  examined  by  looking  at  the  first  few  contributions  from  Stokes 
second-order  wave  theory  where 
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U  _  cosh  *(</■*•:) 
1  iinh  kd 


and 


U,-  4 


3  \wH\  coth  2*W  +  z) 


4  |  L  I  ainh  4  kd 
2w 


(3.6a) 


(3.6b) 


with  d  being  the  water  depth,  k  the  wavenumber  — ,  //the  wave  height,  and  z  the  depth  of  the 
cylinder  below  mean  water  level.  The  two  additional  drag  terms  are 


and 


~  (|cos20|cos0  +  cos 20 1 cos  0|) 
U\ 


2|"7r(  (cos  20 1 cos  20 1) 

*J\ 


(3.7a) 


(3.7b) 


while  the  additional  inertia  term  is 


2|~j  *in  20  (3.7c) 

one  of  which  is  plotted  in  Figure  4.  The  proliferation  of  harmonics  is  quite  clear  and  their 
importance  depends  directly  on  the  ratio  U2/U\.  For  deep  water  waves  of  any  realizable  steep¬ 
ness  this  ratio  is  quite  small.  For  a  wide  range  of  other  wave  conditions  this  ratio  remains  small 
enough  that  the  use  of  Stokes  expansions  can  be  avoided. 


It  is  perhaps  more  interesting  to  examine  the  complementary  problem  of  obtaining  Cm 
and  CD  from  the  same  wave  force  record  using  different  order  Stokes  wave  theories.  The 
integral  least-squares  results,  Eqs.  2.18  and  2.19,  can  be  generalized  to  permit  arbitrary  (but 
orthogonal)  velocities  and  accelerations  as  follows 


_  /„  /mem  v\u\d9 

//*  (u\u\)2d« 

_  L  /**•,»** 

fJw  it2  d9 
•*0 


(3.8) 

(3.9) 
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The  expression  Tor  Cu  is  easier  to  deal  with  and  an  important  point  can  be  made  using  it  alone. 
Upon  substitution  of  u  from  Eq.  3.4  the  expression  for  Cm  becomes 

Cm  -  4  77777  •  C  /«  CL  U,  sin nB)  dB  (3.10) 

in  analogy  with  earlier  expressions.  Next  let  the  non-dimensional  measured  force  be  separated 
into  its  Fourier  components  as  follows 

r 

fm,as(B)  ”  X.  S'n  n 9  +  b"  C0S  «•)•  (3.11) 

Upon  substitution  into  Eq.  3.9  the  result  is 

C* - i  „2  ait/i  +  alVl+aiUi+...  (3.12) 

itlpDlL  I  U„ 

as  one  would  expect.  This  result  is  however  based  on  the  assumption  that  the  harmonics  in  the 
measured  force  are  solely  a  result  of  the  harmonics  in  the  acceleration.  If  harmonics  occur  for 
some  other  reason  then  this  approach  cannot  distinguish  their  origin  and  simply  lumps  the 
effect  into  the  computed  coefficient  and  by  definition  will  give  a  good  (it  to  the  data.  The  prob¬ 
lem  once  again  is  to  attach  any  general  significance  to  the  coefficient  so  obtained.  A  natural 
question  concerns  the  type  of  harmonic  structure  that  can  be  expected  for  various  flows. 

Sarpkaya  largely  neglects  the  higher  harmonics  and  thereby  implies  their  smallness.  On 
the  other  hand  Keulegan  &  Carpenter  devoted  significant  effort  to  categorizing  the  higher  har¬ 
monics  or.  as  they  termed  it,  the  residual  force  The  observed  residual  force  was  greatest  in  the 
neighborhood  of  K  -  IS  which  is  usually  associated  with  asymmetric  vortex-shedding.  The 
peak  residual  force  was  as  much  as  twenty  percent  of  the  peak  force  and  had  components  at 
both  cos  3 B  and  sin  3 B.  It  seems  that  these  variations  in  the  in-line  force  correspond  to  some 
of  the  largest  changes  in  the  transverse  or  lift  force.  Due  to  the  orbital  motion  of  the  fluid,  a 
measurement  of  the  horizontal  component  of  force  /,(#)  on  a  horizontal  cylinder  will  at  times 
contain  a  component  of  the  lift  force  generated  by  the  vertical  component  of  velocity.  Since 
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this  usually  oscillates  at  a  multiple  of  the  wave  frequency,  it  will  appear  as  a  harmonic  in  the 
force  record  and  in  view  of  the  above  equation  it  can  lead  to  an  improper  computation  of  the 
coefficients.  Thus,  higher-order  wave  theories  must  be  applied  with  some  care  in  order  to 
obtain  consistent  results  and  in  order  not  to  mask  some  of  the  physical  differences  between 
various  situations.  By  the  same  token  it  may  be  possible  to  identify  the  various  physical 
processes  by  carefully  analyzing  the  changes  in  harmonic  structure.  It  is  interesting  to  note  that 
Chakrabarti’s  experiments  with  inclined  cylinders  in  waves  produced  the  largest  residual  force 
variations  at  twice  the  wave  frequency  (10,1 1]  as  compared  to  30  for  Keulegan  &  Carpenter.  As 
a  final  comment,  the  sensitivity  to  harmonics  of  single-point  methods  for  obtaining  C»  and  CD 
is  evident  in  the  above  expressions. 

There  are  many  possible  wave-cylinder  configurations,  but  two  particular  cases  are  studied 
here  to  isolate,  as  far  as  possible,  two  different  physical  situations.  The  first  configuration  is  a 
horizontal  cylinder  whose  axis  is  parallel  to  the  wave  crests.  As  in  a  simple  oscillating  flow,  the 
velocities  and  accelerations  are  always  normal  to  the  cylinder  axis  and  it  is  possible  to  select  a 
diameter  that  is  small  compared  to  the  scale  of  the  vertical  gradients  of  velocity  and  accelera¬ 
tion.  Under  these  conditions  one  would  expect  some  correspondence  between  a  wave  flow  and 
the  simple  oscillating  flow.  The  principal  differences  are  due  to  the  eccentricities  of  the  water 
particle  orbits  under  waves,  such  that  the  wake  is  not  always  swept  back  over  the  cylinder.  The 
vertical  cylinder  case  admits  the  axial  variation  of  velocity  and  generally  the  velocities  and 
accelerations  are  no  longer  normal  to  the  cylinder  axis.  On  the  other  hand,  the  wake  is  swept 
back  and  forth  over  the  vertical  cylinder  as  it  is  in  the  one-dimensional  case. 

3.1  Horizontal  Cylinder  —  Even  with  the  available  simplifications  it  is  not  obvious  what 
form  Morison's  equation  should  take  for  wave  flows  about  a  horizontal  cylinder.  This  is  due  to 
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the  second  components  of  velocity  and  acceleration.  It  has  become  customary  to  write 
Morison’s  equation  for  each  component  of  force  using  the  corresponding  components  of  velo¬ 
city  and  acceleration  to  obtain 


and 


/x<«) 


it 

UJ  T 


(3.13) 


where  Um  is  as  yet  an  unspecified  characteristic  velocity  and  the  subscripts  refer  to  the  com¬ 
ponent  directions.  It  is  also  possible  to  write  Morison’s  equation  using  vector  quantities  and 
then  to  take  the  components  of  the  resultant  vector  force  so  that 


A<‘>  -  c°  f  +  T  C“ 


and 


(3.15) 


where 


t  tg)  _  r  .  zL  r  "T 
JzW>  ”  <-o  ni  +  -y  cm  -^7^. 


2ir  Q' 


(3.16) 


7  “  ui  +  wk  (3.17a) 

and  _ 

q  —  VuJ  +  w1,  (3.17b) 

q  -  V(/J+  P  (3.17c) 

The  second  pair  of  equations  for  the  force  components  are  seldom  employed  even  though  they 
appear  somewhat  more  rigorous.  It  must  be  emphasized  that  both  approaches  reduce  to  the 
one-dimensional  case  and  both  give  a  good  fit  to  the  data  but  for  different  values  of  the 
coefficients  Co  and  CM.1  This  is  the  crux  of  the  problem  because  a  common  objective  (and 
one  objective  of  this  investigation)  is  to  determine  what  form  of  the  Morison  equation  is  valid 
for  wave  flows  and  whether  or  not  the  coefficients  obtained  from  a  one-dimensional  experiment 

'The  linearity  of  the  inertia  term  often  result*  in  no  difference  between  the  two  approaches. 
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are  useful  in  a  two-dimensional  wave  flow.  Clearly,  the  choice  between  Eqs.  3.1  and  3.2  or 
Eqs.  3.3  and  3.4  will  have  to  be  made  with  the  aid  of  experimental  results  and  a  careful  compar¬ 
ison  with  one-dimensional  results. 

Although  the  vector  approach  appears  to  yield  the  more  complex  expressions,  the  particu¬ 
lar  case  of  deepwater  waves  produces  a  simpler  expression  for  the  drag  force  and  no  difference 
at  all  between  the  inertia  force  terms.  For  a  horizontal  cylinder  at  a  depth  z  -  —  d(  we  have 
from  linear  wave  theory 

u(t)  =  —jr  e  cos#,  u(i)  —  — e  ‘  sin#  (3.18) 

and 

,  ,  irH  -kd,  .  .  ,  2n2H  -kd, 

*v(t)  »  —jr  e  sm0,  wU)  -  ~~p~  e  cos 9 

TT  H  -kd 

where  S  =»  kx  -  cu/.  By  defining  Um  =  e  '  the  drag  terms  in  Eqs.  3.13  and  3. 14  become 


/to(»)  -  CD  cos0|cos0| 

(3.19) 

and 

/dj(0)  -  CD  sin0|sin0|. 

(3.20) 

The  drag  terms  in  Eqs.  3.3  and  3.4  take  the  alternate  forms 

/dx(0)  “  CD  cos  0 

(3.21) 

and 

foz(9)  -  CD  sin0. 

(3.22) 

The  difference  between  the  two  pairs  of  expressions  is  significant  in  that  both  give  the 
same  peak  drag  force  but  for  the  remainder  of  the  cycle  the  second  pair  predict  larger  drag 
forces  than  the  first.  For  example,  at  9  -  n/4  the  second  approach  gives  a  drag  40  percent 
larger  than  the  usual  approach  for  equal  peak  values.  In  a  situation  where  the  drag  and  inertia 
forces  are  comparable  in  magnitude  this  difference  means  about  a  20  percent  change  in  the 
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predicted  peak  wave  force.  Therefore,  the  proper  selection  of  approach,  i.e.  component  versus 
vector,  has  important  consequences  other  than  the  matter  of  specifying  CM  and  Q>.  An  experi¬ 
ment  that  measures  both  components  of  force  on  a  horizontal  cylinder  under  waves  can  be 
expected  to  resolve  this  question,  but  owing  to  vortex-induced  side  force  this  will  not  be  as 
straightforward  as  it  might  appear  [10,  il]. 

3.2  Vertical  Cylinders  —  As  in  a  planar  flow  the  vortex  wake  of  a  vertical  cylinder  in  waves 
is  swept  back  and  forth  over  the  cylinder.  However,  the  flow  usually  varies  in  direction  with 
time  and  varies  in  magnitude  along  the  span.  The  usual  approach  is  to  write  Morison’s  equa¬ 
tion  in  a  differential  form,  take  the  components  of  velocity  and  acceleration  normal  to  the 
cylinder,  and  then  integrate  over  the  length  of  the  cylinder  to  get  the  total  in-line  force.  The 
key  assumption  is  that  CM  and  CD  can  be  taken  as  uniform  along  the  span  to  get 

/  “  CD  cos  0|cos  0|  J"  ^  a2dz  -  CM  n2/  K  sin  0  J*__  adz  (3.23) 

where 

a  -  cosh  kz  +  tanh  kd  sinh  kz  (3.24) 

for  a  cylinder  extending  from  a  depth  z  ■»  -  zk  to  the  instantaneous  free  surface  rj.  It  is  also 

customary  to  evaluate  the  coefficients  at  the  still  water  level,  which  is  to  say  that  the  charac¬ 
teristic  velocity  U  is  the  horizontal  velocity  at  z  -  0.  The  obvious  difficulty  with  this  approach 
is  its  inability  to  account  for  different  vertical  distributions  which  may  have  the  same  charac¬ 
teristic  velocity  U.  For  example,  two  cylinder-wave  combinations,  one  a  deepwater  wave  and 
one  a  shallow  water  wave,  can  be  selected  to  give  the  same  values  of  K  and  /3  at  z  —  0,  but  the 
vertical  distributions  of  Cm  and  Qj  can  be  quite  different.  The  discrepancy  is  potentially 
greater  when  the  differences  in  orbital  motions  are  considered  as  well.  A  simple  correction  can 
be  obtained  by  assuming  a  linear  distribution  of  CM,  fitting  this  to  the  data  and  then  comparing 
the  linear  distribution  to  a  distribution  constructed  from  the  planar  flow  results.  Since  the  total 
inertia  component  of  force  is  the  same  regardless  of  the  choice  of  description  we  have 
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where 


a  sin  9  dz 


(CWl  +  CM2z)  a  sin  6  dz 


(3.25) 


Cm  -  the  best  fit  constant  coefficient. 

Cmx  —  the  z  -  o  value  of  a  linear  distribution,  and 
•  the  slope  of  the  linear  distribution. 

This  expression  can  be  used  to  estimate  the  slope  of  a  "best  fit"  linear  distribution  by  setting  CM 
equal  to  the  coefficient  obtained  using  Eqs.  3.1  and  3.2  and  setting  CMj  equal  to  the  correspond¬ 
ing  planar  flow  value.  For  simplicity  and  in  view  of  the  approximate  nature  of  this  approach  the 
upper  limit  of  integration  is  taken  as  z  —  o  to  get 

s\  ** 

Cm,  “  <CM  -  CM)  -  (3.26) 

2  1  j  za  dz 

where  the  integrals  are  easily  evaluated  for  each  wave  condition.  A  similar  analysis  for  the  drag 


term  yields 


Cd ,  “  ~  CD  ) 


5° fl2<fe 

J-H 

f  za2dz 


(3.27) 


With  these  expressirns  it  will  be  possible  to  examine,  in  a  rudimentary  way,  the  effects  of  flow 
gradients  along  the  cylinder  axis  upon  the  individual  drag  and  inertia  components  of  the  wave 
force. 


4.0  PRESENTATION  AND  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


The  present  results  were  obtained  for  relatively  low  Reynolds  numbers, 
Re  =  103  -  4  x  103,  and  for  values  of  the  frequency  parameter  in  the  range  /9  —  300-700. 
The  depths  of  the  horizontal  cylinders  varied  between  three  and  twelve  cylinder  diameters 
below  still  water  level  while  the  vertical  cylinders  were  positioned  to  pass  through  the  free  sur¬ 
face  at  all  times  and  terminate  some  distance  (approx.  12  cyl.  dia.)  from  the  bottom  of  the 


22 


■o'. 


tank. 


NRL  MEMORANDUM  REPORT '  420b 

4.1  Horizontal  cylinders.  Despite  the  common  use  of  horizontal  or  near  horizontal  struc¬ 
tural  members  in  many  offshore  platforms,  little  work  has  been  done  regarding  horizontal 
cylinders  under  waves.  From  a  hydrodynamic  viewpoint  the  wave  flow  about  a  horizontal 
cylinder  is  interesting  because  of  certain  similarities  to  the  planar  flow  case.  Namely,  the  flow  is 
uniform  over  the  span  and  is  always  perpendicular  to  the  cylinder  axis.  The  significant 
differences  are  due  to  the  orbital  motions  under  waves  such  that  the  wake  is  not  always  swept 
back  and  forth  over  the  cylinder. 


A  recent  paper  by  Maull  &  Norman  [81  considered  the  horizontal  cylinder  under  waves 
and  represented  the  measured  force  by  RMS  coefficients  defined  as 

C  _  (F0rcev>™' 

'  1/2  PDL(Urms)2 

for  the  force  component  in-line  with  the  wave  direction  and  as 

r’  (Force,.)  rms 

'■  “  1/2  p  DL(Wrmy  2) 

for  the  component  transverse  to  the  wave  direction.  C,  is  often  referred  to  as  a  lift  coefficient 

although  this  is  potentially  confusing  because  the  transverse  force  in  this  case  is  also  of  the 

Morison-type.  The  characteristic  velocities,  U  and  IF,  are  computed  for  the  centerline  depth  of 

the  cylinder. 


Comparisons  between  the  present  results  and  those  of  Maull  &  Norman  are  presented  in 
Figures  5  and  6.  The  coefficients  are  based  on  the  RMS  force  components  at  the  wave  fre¬ 
quency.  The  comparisons  can  only  be  qualitative  since  the  ratio  of  maximum  vertical  velocity 
W  to  maximum  horizontal  velocity  V  in  the  present  experiments  ranges  from  W/U  —  0.9  to 
W/U  -  1.09  while  Maull  &  Norman’s  data  were  obtained  for  W/U  -  0.81-0.84.  Also,  their 
velocities  and  accelerations  were  computed  from  Stokes  second  order  theory  whereas  the 
present  values  were  obtained  from  the  linear  theory.  Maull  &  Norman  do  not  give  a  range  of /3 
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values  although  the  experimental  description  suggests  that  they  should  be  close  to  the  present 
values  of  0  —  4S0  and  /3  —  630.  In  spite  of  these  differences  the  results  of  the  two  sets  of 
experiments  agree  quite  well.  Also  included  in  the  figures  ape  Maull  &  Norman’s  representa¬ 
tions  of  the  predicted  RMS  coefficients  based  on  Sarpkaya’s  planar  flow  results.  A  substantial 
overprediction  is  observed  although  it  must  be  noted  that  as  W /  U  decreased  Maull  &  Norman 
found  that  the  RMS  wave  forces  tended  toward  the  one-dimensional  results  as  might  be 
expected. 

The  total  in-line  and  transverse  RMS  coefficients  for  the  present  results  are  compared  in 
Figure  7  to  each  other  and  to  the  planar  flow  results  of  Sarpkaya  [5]  and  of  Bearman,  Graham 
&  Singh  [6],  The  discrepancies  between  the  two  sets  of  one-dimensional  results  are  unex¬ 
plained  but  again  either  set  will  significantly  over-predict  the  RMS  wave  force.  It  should  also 
be  noted  that  the  RMS  transverse  or  lift  force  is  consistently  less  than  the  RMS  in-line  force 
for  equivalent  values  of  the  Keulegan-Carpenter  number. 

The  individual  coefficients  Of  and  CD  can  be  determined  once  a  suitable  form  of 
Morison’s  equation  is  adopted.  The  choices  were  described  in  Section  3.1  and  amounted  to 
either  a  scalar  application  of  Morison’s  equation  to  the  components  of  velocity  and  acceleration 
or  a  vector  approach  using  the  vectors  of  velocity  and  acceleration.  The  best  fit  of  each  method 
to  a  typical  measured  force  record  is  presented  in  Figure  8  where  it  is  evident  that  the  vector 
approach  is  indeed  more  appropriate  and  suggests  the  use  of  Eqs.  3.3  and  3.4  for  horizontal 
cylinder.  In  every  case  Eqs.  3.1  and  3.2  produce  drag  coefficients  that  are  smaller  (approx.  5- 
10%)  than  the  drag  coefficients  obtained  from  Eqs.  3.3  and  3.4.  The  inertia  coefficient  is,  of 
course,  the  same  in  each  case.  Another  significant  feature  of  the  vector  approach  is  the  linear¬ 
ity  of  the  resulting  force  components  in  deep  water  which  will  permit  stochastic  analysis  of  hor¬ 
izontal  members  without  artificial  linearization. 
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The  least-squares,  best  fit  inertia  coefficients  are  presented  in  Figure  9  along  with  the 
one-dimensional  results  of  Sarpkaya  and  of  Bearman  et  al.  As  might  be  anticipated  from  the 
RMS  coefficients,  the  wave-induced  intertia  forces  are  less  than  either  set  of  planar  flow  results 
but  appear  to  follow  the  same  trend  of  decreasing  CM  with  increasing  K.  On  the  other  hand  the 
results  for  the  drag  coefficient  exhibited  a  marked  departure  from  the  one-dimensional  results 
as  shown  in  Figure  10.  Not  only  are  the  drag  coefficients  less  than  their  planar  flow  analogues 
but  they  also  do  not  contain  the  "hump"  which  is  characteristic  of  CD  in  this  range  of  K.  There 
is  a  suggestion  that  the  present  results  might  eventually  merge  with  the  planar  flow  results  of 
Bearman  et  al.  at  larger  values  of  K  and  this  ought  to  be  investigated  further. 

The  difference  in  behavior  between  the  wave  drag  coefficients  and  the  planar  flow  drag 
coefficients  must  be  explained  on  the  basis  of  the  orbital  motions  of  the  fluid.  Physically,  this 
is  satisfying  because  the  present  range  of  Keulegan-Carpenter  number  corresponds  to  a  rela¬ 
tively  small  number  of  vortices  being  shed  during  a  flow  cycle.  If  these  few  vortices  are  rectil- 
inearly  swept  back  over  the  cylinder  then  their  influence  on  the  nascent  vortices,  thence  the 
force,  is  likely  to  be  much  larger  or  at  least  much  different  than  a  few  vortices  which  are  swept 
away  during  the  orbital  motion  under  waves.  In  a  sense,  the  drag  on  a  horizontal  cylinder 
under  waves  is  more  akin  to  the  steady  flow  drag  than  to  the  1-D  unsteady  drag.  Thus,  the 
"hump"  behavior  will  not  occur  for  horizontal  cylinders  and  the  drag  coefficient  appears  to 
monotonically  approach  the  steady-state  drag  coefficient  as  the  Keulegan-Carpenter  number 
increases.  This  limit  has  been  demonstrated  by  Bearman  et  al.  for  a  variety  of  body  shapes  in 
planar  flows  at  large  K  [6]. 

4.2  Vertical  cylinders— The  present  results  are  listed  together  with  corresponding  one¬ 
dimensional  data  in  Table  1.  The  one-dimensional  values  were  obtained  from  Sarpkaya  [S] 
because  his  data  are  quite  extensive.  Even  so,  it  was  often  necessary  to  interpolate  and,  in 
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Table  1  —  Vertical  Cylinder  Results 


*>,-0 

0>,-o 

Total 

Meas. 

RMS 

Force 

Total 

Pred. 

RMS 

Force 

RMS 

Inertia 

Force 

Ratio 

RMS 

Drag 

Force 

Ratio 

Q, 

Co 

w/u 

du  1 

*Lo 

cL 

5.8 

701 

1.16 

2.18 

0.58 

0.05 

0.76 

0.11 

0.94-1.0 

0.81 

0.08 

8.2 

454 

1.44 

2.58 

0.55 

0.18 

0.66 

0.36 

0.62-0  98 

0.31 

0.16 

11.3 

573 

0.74 

1.36 

1.01 

0.12 

0.65 

0.26 

0.81-1.0 

0.87 

0.06 

11.4 

619 

0.72 

1.48 

I.Ot 

0.14 

0.72 

0.30 

0.87-1.0 

1.11 

0.06 

12.7 

395 

1.12 

2.56 

0.92 

0.17 

0.65 

0.29 

0.52-0.95 

0.31 

0.17 

13.6 

526 

0.70 

1.80 

0.86 

0.16 

0.58 

0.36 

0.74-1.0 

0.81 

0.08 

13.7 

451 

1.06 

2.42 

0.93 

0.26 

0.62 

0.63 

0.62-0.98 

0.51 

0.18 

18.0 

312 

1.16 

3.26 

0.90 

0.22 

0.68 

0.43 

0.41-0.87 

0.21 

035 

some  instances,  to  extrapolate  the  one-dimensional  data  so  that  the  comparisons  should  only  be 
viewed  as  qualitative.  Had  the  one-dimensional  results  of  Bearman  et  al.  [6]  been  used,  the 
comparisons  between  wave  results  and  planar  flow  data  would  be  more  favorable  but  the  same 
general  discrepancies  would  occur.  The  total  wave  force,  as  measured  by  the  RMS  value,  is 
again  over-predicted  by  the  planar  flow  results  but  it  is  not  obvious  whether  this  is  due  to  inap¬ 
propriate  coefficients  of  drag  and  inertia  due  to  the  assumption  of  uniform  coefficients  over  the 
span  or  due  to  the  combination.  To  begin  to  resolve  this  question  the  individual  RMS  com¬ 
ponents  of  drag  and  inertia  were  compared  to  their  predicted  values  and  the  ratios  are  listed  in 
Table  1.  Evidently,  the  overprediction  is  principally  in  the  drag  contribution  since  the  ratio  of 
measured  and  predicted  RMS  inertia  forces  is  close  to  one.  In  fact,  the  values  are  close  enough 
to  suggest  that  a  correction  for  the  distribution  of  CM  may  be  able  to  account  for  the 
discrepancy.  A  simple  estimate  of  the  correction  can  be  obtained  by  assuming  a  linear  distribu¬ 
tion  of  Ca 4,  fitting  this  to  the  data  and  then  comparing  the  linear  distribution  to  a  distribution 
constructed  from  the  planar  flow  results  as  discussed  in  Section  3.2.  Plots  which  are  typical  of 
the  results  that  can  be  obtained  are  presented  in  Figure  11.  Since  the  inertia  force  decays 
essentially  exponentially'  with  depth  in  the  present  experiments  it  is  clear  that  the  actual  CM  dis¬ 
tribution  is  well  represented  by  a  distribution  constructed  from  one-dimensional  flow  results. 
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Alternatively,  a  uniform  value  can  be  assumed  and  the  error  in  doing  so  can  be  estimated  from 
one-dimensional  results  using  Eq.  3.26.  A  plot  of  this  estimate  of  the  error  is  given  in  Figure 
12. 


A  similar  analysis  of  the  drag  component  is  far  less  fruitful  because  of  the  large  disparities 
between  the  observed  and  the  predicted  drag.  The  disparity  in  each  case  can  be  attributed  to 
either  the  variation  in  the  velocity  along  the  span  or  to  the  variation  in  the  flow  direction  during 
the  cycle  or  to  both  of  these.  During  the  present  experiments  the  slopes  of  the  velocity  varia¬ 
tion  and  the  eccentricities  of  the  water  particle  orbits  varied  from  test  to  test.  Unfortunately 
they  often  varied  together  and  so  it  is  not  possible  to  separate  the  two  effects.  The  eccentrici¬ 
ties  W/U  and  velocity  variations,  characterized  by  the  slope  at  z  *  0,  are  listed  in  Table  1.  The 
apparent  trend  is  that  as  the  velocity  variation  decreases  and/or  the  orbits  flatten  the  drag  tends 
toward  the  one-dimensional  result.  With  further  study,  or  possibly  further  analysis  of  existing 

/ 

wave  force  data,  this  approach  could  provide  a  basis  for  modifying  one-dimensional  results  for 
use  in  wave  force  prediction.  It  would  also  be  possible  to  establish  error  bounds  for  the 
assumption  of  a  constant  Cp. 

The  final  entry  in  Table  1  is  the  measured  RMS  lift  coefficient,  which  in  this  case  is 
purely  a  vortex  effect  and  is  not  described  by  Morison's  equation.  In  the  same  sense  that  the 
velocity  variation  and  orbit  eccentricity  influence  the  drag,  one  intuitively  expects  that  the 
vortex-induced  lift  will  be  similarly  influenced  and  indeed  this  appears  to  be  the  case.  An 
important  event  which  was  observed  during  the  experiments  but  masked  by  the  RMS 
coefficient  of  lift  was  a  significant  peak  in  the  lift  that  occurred  near  the  crest  of  the  wave.  The 
peak  resultant  force  in  those  cases  was  much  larger  than  expected  from  any  combination  based 
on  the  RMS  values  of  the  in-line  and  the  lift  force.  This  has  been  discussed  elsewhere  (5,  7) 
and  is  only  mentioned  here  as  a  reminder. 


27 


RAMBERG  AND  NIED2WECKI 


5.0  SUMMARY  AND  CONCLUSIONS 

The  results  of  the  present  investigaton  have  shown  that  the  resolution  of  a  Morison-type 
approach  to  wave  force  prediction  can  be  improved  if  consistent  methods  of  computation  are 
employed  along  with  some  accounting  for  the  major  differences  between  wave-induced  vortex 
flows. 


In  the  case  of  a  horizontal  cylinder  in  waves  the  horizontal  and  vertical  components  of  the 
wave  force  must  be  obtained  from  a  modified  form  of  Morison’s  original  equation  whereby 
each  component  is  derived  from  an  equation  written  for  the  instantaneous  total  force.  The 
modified  equations  are  linear  and  are  therefore  readily  employed  in  stochastic  analyses.  The 
force  transfer  coefficients  CM  and  CD  for  the  horizontal  cylinder  cannot  be  obtained  from  one¬ 
dimensional  planar  flow  results  except  possibly  at  large  Keulegan-Carpenter  numbers.  The 
planar  flow  coefficients  can  overpredict  the  actual  wave  force  on  a  horizontal  cylinder  in  waves 
by  as  much  as  100  percent.  The  present  results  along  with  some  simple  physical  arguments 
suggest  that  the  drag  coefficient  on  the  horizontal  cylinder  rapidly  approaches  the  steady  flow 
value  for  increasing  Keulegan-Carpenter  number  and  does  not  exhibit  the  "hump’  behavior 
which  is  characteristic  of  unsteady  planar  flows. 

The  wave  forces  on  vertical  cylinders  were  decomposed  into  inertial  and  drag  contribu¬ 
tions  for  which  vertical  distributions  were  estimated.  It  was  shown  that  the  usual  assumption  of 
uniform  CM  and  CD  can  lead  to  significant  errors  but  that  a  distribution  of  CM  based  on  one¬ 
dimensional  results  adequately  accounts  for  the  wave-induced  inertia  force.  On  the  other  hand, 
the  one-dimensional  flow  drag  coefficients  greatly  overpredict  the  wave-induced  drag  and  the 
amount  is  attributed  to  the  "steepness”  of  the  velocity  gradient  along  the  span  and  to  the  eccen¬ 
tricity  of  the  water  particle  orbits. 
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The  implications  of  the  present  findings  for  the  wave-induced  forces  on  inclined  cylinders 
and  for  the  forces  on  cylinders  in  short-crested  seas  are  straightforward.  When  the  various 
hydrodynamic  complexities  associated  with  wave  flows  were  isolated  as  much  as  possible,  it  was 
found  that  a  reduction  in  force  occurred  as  compared  to  one-dimensional  results.  Thus  it  may 
be  expected  that  situations  where  all  or  most  complexities  occur  together,  such  as  inclined 
cylinders  or  short-crested  seas,  a  further  reduction  in  wave  force  is  likely.  Predictions  based  on 
results  from  other  configurations,  particularly  one-dimensional  flows,  will  be  very  conservative. 
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APPENDIX  A 

The  Least-Squares  Error  Method 

There  are  a  number  of  methods  that  can  be  used  to  determine  the  coefficients  Cw  and  Cp 
for  a  particular  form  of  Morison's  equation  and  a  certain  set  of  data.  The  most  common  and 
widely  accepted  methods  involve  Fourier  analysis  or  curve-fitting  or  both.  A  complete  descrip¬ 
tion  of  Fourier  analysis  can  be  found  in  most  standard  math  texts  but  descriptions  of 
curvefitting  methods  are  usually  restricted  to  the  simplest  case(s).  For  this  reason  a  brief 
review  of  the  general  least-squares  curvefit  is  given  here  and  the  reader  is  directed  to  a  book  by 
Wolberg  for  a  full  description  [17], 

The  first  step  in  any  curvefit  is  to  assume  that  a  "true"  relationship  exists  in  the  form 

V  “/(([.  (i.  atj  ...am)  (Al) 

where 


n  is  the  dependent  variable. 

( ,  are  the  n  true  independent  variables, 
a ,  are  the  m  true  coefficients,  and 
/(  )  is  the  true  relationship. 

The  objective  is  to  construct  the  best  representation  of  this  relationship  using  measured  values 
of  the  independent  variables  Xk,  and  measured  values  of  the  dependent  variable  K*.  The  con¬ 
structed  equation  is  written  as  follows 

vk  -  /<*,,.  xlt,  xh  . . .  xm,  a,.  a2.  a  t.  am)  (A2) 

yk  is  the  k'h  computed  value  of  the  dependent  variable. 
xit  is  the  k'h  computed  value  of  the  i'"  independent  variable,  and 
a„  are  the  computed  coefficients. 

The  "best  constructed  equation"  is  defined  as  the  equation  whose  values  for  the  coefficients  am 
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satisfy  some  minimum  error  criterion  with  a  particular  set  of  observed  values  Yk  and  Xkl.  It  is 
customary  to  minimize  the  differences  between  the  observed  values  and  the  calculated  values  in 
the  least  square  sense.  To  do  this  we  first  must  define  the  residuals, 


—  ^k  1 

0  —  y  „  I  —  observed  value  -  calculated  value  (A3) 

Kxkl  —  Ak,  ~  Xkl  J 

and  the  weights. 


wyk  ~  1/** 
w,  =  l/tr  1 

*ki  *ki 


(standard  deviation)  1 


(A4) 


The  latter  statistical  weighting  is  adopted  so  that  the  residuals  can  be  measured  in  units  of  stan¬ 
dard  deviation  or,  in  other  words,  units  of  equal  uncertainty.  The  computed  values  of  a„  are 
sought  which  minimize  the  sum  S  of  the  squares  of  the  weighted  residuals 


5 


I 

ft-  1 


R*  +  T  h\  Rx2 

I 


(AS) 


The  minimization  of  5  is  the  general  least  squares  curvefit  The  solution  procedure  depends 
greatly  on  the  form  of  the  "true"  relationship  /and  ultimately  the  form  of  S. 


A  few  points  are  worth  noting.  First,  the  least  squares  curvefit  does  not  determine  if  the 
assumed  relationship  is  true  or  even  if  it  is  a  good  representation  The  "goodness"  of  the 
assumed  function  must  be  judged  by  some  other  means  such  as  its  ability  to  match  all  or  cer¬ 
tain  of  the  observed  values  of  the  dependent  variable  or  perhaps  its  ability  to  produce  certain 
values  of  the  coefficients.  The  question  of  "goodness"  becomes  the  question  of  the  validity  of 
Morison's  equation  for  describing  wave  forces.  In  the  case  of  wave  forces,  the  judgement  of 
"goodness"  of  Morison's  equation  can  be  made  very  difficult  by  the  fact  that  the  wave  force  is 
inherently  periodic  and  the  equation  itself  can  be  viewed  as  the  first  few  terms  of  a  Fourier 
expansion  Thus,  a  reasonably  good  fit  can  always  be  obtained  quite  apart  from  any  physical 
considerations  of  the  flow 
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The  second  comment  concerns  the  contribution  of  the  independent  variable  residuals  to 
the  sum  of  the  weighted,  squared  residuals.  In  most  experiments  the  uncertainties  in  the 
independent  variables  are  much  smaller  than  the  uncertainty  in  the  dependent  variable,  i.e. 
**>  «  w,,  ,  so  that  the  sum  S  reduces  to 

Jk 

s  -  £,  wytRy]  (A6) 

*-l 

and  if  wyk  —  constant  =  w  then 

5  -  w  f  ( Kobserved  -  Kcllc)2  =wE  (A7) 

*-i 

Now  if  the  parameters  a„  appear  linearly  in  the  assumed  function,  /,  then  the  least  square  solu¬ 
tion  is  obtained  from  the  equations  resulting  from 

-  0.  (A8) 

da„ 

This  is  the  commonest  least-square  solution  and  the  one  presented  in  most  texts.  The  point  to 
be  made  here  is  that  this  simple  curvefit  is  often  applied  to  the  wave  force  problem  without 
regard  to  the  implications.  Morison’s  equation  relates  the  wave  forces  to  the  wave-induced 
velocities  and  accelerations  but  these  quantities  are  rarely  measured.  Instead  the  velocities  and 
accelerations  are  computed  from  a  wave  theory  and  the  observed  surface  profile.  This  can  be 
viewed  as  the  assumption  of  another  ’true'  relationship  or  as  the  introduction  of  significant 
residuals  in  the  independent  variables.  In  the  first  instance  the  'goodness'  of  Morison's  equa¬ 
tion  cannot  be  separated  from  the  "goodness"  of  the  wave  theory.  In  the  second  instance  the 
simple  curvefit  procedure  no  longer  applies  and  the  general  form  of  S  must  be  employed  with 
the  introduction  of  more  parameters.  The  general  approach  will  not  yield  new  or  even  very 
different  answers  but  it  would  provide  a  rigorous  basis  for  assessing  the  errors  and  ultimately 
the  degree  of  confidence  in  the  results  obtained. 

The  final  comment  concerns  the  use  of  an  integral  representation  of  the  error  function  to 
be  minimized.  These  expressions  are  generally  of  the  form 
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and  the  least-square  solution  is  again  found  by 


provided  that  E  is  linear  in  the  unknown  coefficients  am.  Obviously,  this  error  criterion  is  the 
limiting  case  of  the  earlier  discrete  summation  and  is  strictly  defined  only  when  W(t)  and  R(t) 
are  continuous  or  piece-wise  continuous  functions.  If  one  wishes  to  apply  the  integral  error  cri¬ 
terion  to  a  discrete  function  then  some  care  must  be  taken  in  evaluating  the  integral  so  that 
arbitrary  statistical  weightings  are  not  introduced.  For  example,  the  integral  may  be  evaluated 
by  simple  rectangular  segments  to  obtain 


E  -  /  r  w(t)  R,*(t)  dt  =  £  fVUt)  R2(0  At ’• 


or  by  Simpsons  rule  to  obtain 


E  =  j  Arimr,)  R2X  +  WU„)  R2)  +  y  M  'l  w(r2,)  •  R l 


+  ~  M  I  w(r2l_,)  Rl_ ,  (A12) 

J  f-2 

which  on  comparison  with  the  usual  discrete  form 

-  £  wF  R,2  (A13) 

»— i 

shows  that  the  choice  of  method  to  evaluate  the  integral  can  produce  very  different  statistical 
weightings  which  will  lead  to  different  least-squares  solutions.  The  integral  operator  has  the 
advantage  of  simplifying  many  calculations  so  it  is  of  interest  to  know  the  conditions  that  lead 
to  the  proper  discrete,  least-squares  solution.  Comparing  Eqs.  All  and  A13  and  keeping  in 
mind  that  the  solution  will  be  given  by  Eq.  A10,  it  is  clear  that  uniform  statistical  weighting  and 
uniform  sampling  intervals  A  twill  permit  the  interchangable  use  of  integral  and  discrete  expres¬ 
sions  provided  that  the  integral  is  evaluated  simply  (i.e.  —  by  means  of  a  rectangular  or  a  tra¬ 
pezoidal  rule). 


APPENDIX  B 

Experimental  Systems  and  Methods 


The  experiments  were  performed  in  a  100  ft  (30.5m)  long,  4  ft  (1.2m)  wide,  and  6  ft 
(1.8m)  high  wave  tank  at  the  Naval  Research  Laboratory  [18]  (see  Figure  Bl).  The  tank  was 
filled  to  a  3-1/2  ft  (1.1m)  depth  during  the  tests.  A  one  inch  (2.54cm)  diameter  cylinder,  3  ft 
(0.9m)  long,  was  placed  in  a  U-shaped  frame  with  two-component  force  gages  at  each  end.  The 
frame  could  be  positioned  so  that  the  cylinder  was  either  vertical  or  horizontal.  In  the  horizon¬ 
tal  configuration  the  bulk  of  the  support  frame  was  above  water  and  the  entire  assembly  was 
raised  or  lowered  to  achieve  the  desired  cylinder  depth.  For  the  vertical  cylinder  tests  the  U- 
shaped  frame  was  fastened  to  a  vertical  strut  near  one  side  of  the  channel  so  that  the  measure¬ 
ment  cylinder  was  centered  in  the  flume.  The  gap  between  the  measurement  cylinder  and  the 
nearest  vertical  support  member  was  greater  than  12  cylinder  diameters.  End  effects  1 14]  were 
not  examined  in  this  study  but  will  be  the  subject  of  subsequent  work  on  horizontal  cylinders  in 
waves. 

The  force  gages  were  assembled  as  square-section  cantilever  beams  with  full-bridge 
straingage  circuits  arranged  to  sense  orthogonal  components  ot  only  the  force  applied  to  normal 
to  the  beam's  end.  The  cylinder,  force  gage,  and  U-frame  assembly  was  statically  calibrated  by 
applying  known  deadweights  at  various  positions  and  directions  along  the  cylinder.  The  linear¬ 
ity  of  the  output  as  well  as  the  proper  partitioning  of  the  load  between  the  sensors  at  each  end 
is  indicated  in  Figure  B2.  From  these  tests  the  equivalent  load  for  a  precision  shunt  resistor 
across  one  arm  of  the  bridge  was  also  determined  so  that  the  resistor  could  be  switched  into  the 
circuit  for  later  check  calibrations.  The  signals  from  the  four  DC  bridge  circuits  were  amplified 
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and  filtered  (Preston  DX-I  amplifiers)  before  plotting  on  a  strip  chart  recorder.  The  dynamic 
response  of  the  force  measurement  and  recording  system  was  checked  by  pendulum  oscillations 
of  the  deadweights.  The  phase  shifts  between  the  applied  signals  and  the  recorded  results  were 
undetectable  (less  than  0.01T  for  0  <  1/T  <  2hz  and  less  than  0.05T  for  2  <  1/T  <  lOhz). 


After  recording  the  data  on  strip  charts  several  consecutive  cycles  were  digitized,  by  hand, 
and  stored  on  tape  for  processing  with  a  desktop  computer  (Tektronix  40S1). 


The  period  and  height  of  the  incident  wave  were  set  by  the  motions  of  a  mechanical, 
bulkhead-type  wave  generator  located  at  one  end  of  the  channel  [18).  For  the  present  range  of 
conditions,  98  percent  or  more  of  the  generated  wave  energy  was  absorbed  by  a  porous,  sloping 
beach  at  the  other  end  of  the  channel.  The  instantaneous  water  elevation  at  the  cylinder  was 
recorded  and  digitized  in  the  same  manner  as  the  forces  with  a  capacitance-type  wave  gage  as 
the  sensor.  The  measurement  system  was  calibrated  before  each  run  and  checked  after  each 
run.  The  accuracy  of  all  measurements  is  believed  to  be  better  than  ±2  percent. 


The  force  transfer  coefficients  were  obtained  from  the  digitized  wave  force  record  by  the 
least-squares  criterion, 

S-£  (fma<  -  fcaU  >2  (Bl) 

1—1 

subject  to 


as  _  _  o 

acm  bcd 

where  N  is  the  number  of  data  points  in  the  cycle  which  varied  from  20  to  35. 
induced  velocities  and  accelerations  were  computed  using  linear  wave  theory. 


(B2) 
The  wave- 


37 


C* 


■% 


RAMBERG  AND  NIEDZWECKI 


I 


Harmonic  comribuiions  10  th 
2nd  order  Stokes  expansion 
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Keu  1  egarv-Carpenter  Number,  K 

Figure  6  -*  The  transverse  RMS  wive  force  coefficients  (evaluated  for  the  fundamental 
frequency)  for  a  horizontal  cylinder  in  waves.  Present  data:  IS  ~  430.  fl  “  610.  •. 
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Figure  8  -  Two  leasi-squares  besi  fils  lo  a  measured  wave  force  on  a  horizontal  cylinder 
Usual  Morison  approach.  —  (Eqs  3  13  &  3  14).  vector  approach.  -  -  •  (Eq  3  15  &  3  16) 


Coef f 1 c 1 ent 


Keulegin-Carpenter  Number,  K 

Figure  10  —  The  in-line  end  transverse  drag  coefficients  for  horizontal  cylinders  in  waves  and 
for  corresponding  one-dimensional  oscillatory  flows  Present  data,  in-line.  O;  transverse.  •. 


Figure  1 1  —  A  comparison  between  a  best-fit"  linear  distribution  of  C  n  for  a  vertical  cylinder  in  waves  and 
a  distribution  constructed  from  one-dimensional  oscillatory  (low  results;  see  Eq.  (3  25). 
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f  the  error  in  the  computed  inertia  force  on  a  vertical  cylinder  in 
constant,  one-dimensional  flow  values  of  Cw  are  employed. 


A 


Figure  B1  —  A  photograph  of  the  wave  lank 
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